ON THE CAUCHY PROBLEM FOR A BOUSSINESQ TYPE 

SYSTEM 



F. LINARES, M. PANTHEE, AND J. DRUMOND SILVA 

Abstract. We consider the initial value problem (IVP) associated to a Boussinesq 
type system. After rewriting the system in an equivalent form of coupled KdV-type 
equations, we prove that this is locally well-posed in (iJ s (R 2 )) 4 , s > 3/2, using sharp 
smoothing estimates. Consequently we obtain the local well-posedness result for the 

' 

original system in H s x V s+ for s > 3/2 (see below for the definition of V s ). 



Oh 

< 



1. Introduction 
We consider the following Boussinesq type system 

'ry t + A$-fA 2 $ = -eV- [r)((d Xl <S>y,(d X2 <S>y)], 

{t, X) — [t, X\,X2) £ 

® t + v -fi(a- 1/2) A V = --^((d Xl $y +1 + (d X2 <s>y +1 ), 



)l+2 



(i.i) 

where e is the amplitude parameter, fx is the long-wave parameter and a is the Bond 
number. This system was derived in [S] as a model for the evolution of three dimensional 
; long water waves with small amplitude, where $ represents the non dimensional velocity 

^ ! potential at the bottom z = and rj represents the free surface elevation. 

Well-posedness of the associated Cauchy problem and stability of solitary waves 
for this system was studied in [9]. In this paper the author proved that the initial 
value problem (IVP) associated to (11.11) is locally well-posed for given initial data in 
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H S (M. 2 ) x V S+1 (M 2 ), s > 2, where H S (R 2 ) is the usual L 2 -based Sobolev space and 
V s (IR 2 ) is the Hilbert space defined by 



V S (R 2 ) := {/ G S'(K 2 ) : V^Af G H 3 - 1 }, 



with the corresponding norm ||/||v s (r 2 ) = \\\ / —Af\\H 3 - 1 (R 2 )- The idea used in [H] is 
based on a change of variables by taking u = d Xl Q and v = d X2 $>, thus enlarging the 
system to a new one for U = (j],u,v), and then employing Kato's approach for the 
generalized Korteweg-de Vries (KdV) equation. The author also proves the existence 
of global weak solutions in time, for small initial data, by considering a variational 
approach, which involves the Hamiltonian defined in the energy space H 1 x V 2 . 

In this work we aim to improve the local results obtained in [9]. For the sake of 
simplicity we consider yU = 6, a = 2/3 and e = 1 in (11. ip and we start by diagonalizing 
the constant coefficient linear part of the resulting system 



T) t + A$ - A 2 $ = -V ■ [r){® p xi , 
$ t + r? _ Ar? = __l_ ( $ S +i + $P+i) 

to obtain the following equivalent system 



;i.2) 



u t + iv=A(l - A)u = -JVi(u,v) - N 2 (u,v), 



;i.3) 



v t - i v /z A(l - A)v = JVi(u, f ) - iV 2 (w, f ), 
where the nonlinearities are given by 

N x {u, v) = ^{<9 2l [{V=£(u - v)}(u Xl + v xi y] + d X2 [{V^A(u - v)}(u X2 + v X2 Y]}, 

N*(u,v) = 2$+r ) [(u xl + v xi y+ 1 + (u X2 + v X2 y +1 ]. 

(1.4) 

In the diagonalization process the relation between the original functions (rj, $) and 
the new functions (u, v) is given by 

2 ^ +2 ^' (1.5) 
so that (rj, $) can be recovered from (u, v) by 



?7 = i^—Au — i\/—Av, 

(1.6) 
$ = u + v. 
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This change of variables transforms the original system ( II. 2p . with different differ- 
entiability requirements for (r),<&), into a more symmetric form (II. 3p of two coupled 
KdV-type equations with the same third order dispersive operator in the spatial vari- 
ables i\J — A(l — A) for (u, v). The main difficulty now stems from the nonlinearities 
(II. 4p . involving powers of the derivatives of the unknowns. To deal with this issue, 
we differentiate the diagonalized system in x\ and 22, enlarging it to a system of four 
equations for the partial derivatives of the unknowns (u Xl ,u X2 ,v xl ,v X2 ) (see equation 
(13. ip below) thus turning the nonlinearities into standard KdV-type ones. It is this 
final system that we then solve, using classical methods (see |6j) involving Strichartz, 
maximal and smoothing estimates, which are developed in the next section of this 
paper. 

In the third section we then prove the local well-posedness result that constitutes 
the main theorem of this paper. 

Theorem 1.1. Let p = 1 and s > 3/2, then for given data (770, $0) £ H S (K 2 ) x 
V S+1 (IR 2 ), there exist a time T = T(||r7 ||_H" s (R 2 ) + II^oIIv+^r 2 )) > 0, a space Vf, contin- 
uously embedded in C([0, T] : H S (M. 2 ) x V S+1 (IR 2 )) and a unique solution (77, $) G "ty s T to 
the IVP associated to (II. ip . 

Throughout the rest of this paper, we consider p = 1 in ( 11.41) . 
Remarks: 

i) Although we only treat the case p = 1, the same method can be adapted for any 
integer value of p > 1. The regularity for which LWP can be achieved will naturally 
depend on p. 

ii) This result not only extends the LWP claim in [9] but it also provides a complete 
and detailed proof of the equivalence between the original system and the one 
obtained through differentiation, which seems to be missing there. 

Notation: Finally we introduce some notation that we will use in this work. We use 
/ to denote the Fourier transform of / and H s = (1 — A)~ S / 2 L 2 to denote the L 2 -based 
Sobolev space of order s. The Riesz potential of order — s is denoted by D s x = (—A) 5 / 2 , 
so that one has D x = y/— A with Fourier symbol given by |£|. We also define the Riesz 
transforms Ri, I — 1, 2 via the Fourier transform by Ri(f) = ( — ifnf) ■ 
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For / : W 1 x [0, T] -}Kwe define the mixed L^L^-norm by 




(/ T iL lf{x,t)lPdx ) 



dt 



) 



with usual modifications when p = oo. An analogous definition is used for the other 
mixed norms L^L^, with the order of integration in time and space interchanged. We 
replace T by t if [0,T] is the whole real line IR. 

We use c to denote various constants whose exact values may vary from one line to 
the next. We use A < B to denote an estimate of the form A < cB for some c, and 
A ~ B if A < B and B < A. Also, we use the notation a+ to denote a + e for e > 0. 



In this section we present some estimates satisfied by the group associated to the 
linear part of the diagonalized system f ll.3p . These estimates are very similar to the 
analogous ones in [7J, and are actually simpler here due to the better symbol of the 
phase of the linear propagator. 

We thus consider the following pair of linear equations 



whose solution is described by the unitary groups 

is the Fourier symbol given by </>(£) = |£| 3 + |£|. Note that, as in [7], this symbol 
is a radial function, whose specific Fourier transform properties are a crucial ingredient 
exploited in the proofs that follow. We will only derive estimates associated to the 
unitary group U := U + , as the estimates for U~ are the same. 
We start with the smoothing estimate. 

Theorem 2.1. Let T > and {Q a } ae z 2 oe a family of non-overlapping cubes of unit 
size so that IR 2 = Li a ez 2 Qa- Then the following estimate holds 



2. Linear Estimates 



w t ±i^K{\ - A)w = 
w(x, 0) = w (x), 



(2.1) 




(2.2) 



where the implicit constant does not depend on T. 
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Proof. Without loss of generality, we can consider the cubes Q a = Q := {x : \x\ < 1} . 
We have that the gradient of the Fourier symbol |V0(£)| = 3|£| 2 + 1 is always positive. 
Therefore from Theorem 4.1 in [3], it follows that 

\\DlU{t)w a \\ L%M < ( / ^J^I^(OI 2 <T S IkolUl- (2.3) 

□ 

Next we present the maximal function estimate. In the following theorem, the Fourier 
symbol will be 0(0 = |£| 3 + |£|, with f e M. n . Even though we only need the 2- 
dimensional case, this result can be established n dimensions. 



Theorem 2.2. Let {Q a }aez n be the mesh of dyadic cubes in M. n . Then for any s > A 



3n 

in uj uyuuiu cuuts in jr*. . ± ncn jui any s >» 

and T > 0, the following estimate holds 



( V sup sup \U(t)w (x)\ 2 Y 2 < (l + T^)\\w \\ H s. (2.4) 

V a&n \t\<Tx£Q a / 

The proof of this theorem is totally analogous to the corresponding one in [7J, relying 
heavily on the particular representation of the Fourier transform of radial functions in 
terms of Bessel functions. So we omit this proof here, referring the reader to the fully 
detailed proof in that paper, observing only that we now need to consider the phase 
0(0 = I0 3 + 10 (instead of 0(0 = e|0 3 - 10 there). 

As this representation of the Fourier transform of a radial function will also be 
required in the remaining part of this section, we gather some properties satisfied by 
it. The Fourier transform of a radial function /(|ar|) = f(s) is still a radial function 
and is given by (see [10] ) 

POO 

f{r) = fm) = r~— / f(s)J^(rs)s*ds, (2.5) 

Jo 

where J m is the Bessel function defined by 

(r/2) m f 1 i 1 

Mr) = = F LL /,„, x / e™{l-s 2 T-*ds, for m > --. (2.6) 
1 [m + 1/2J7T2 J -i * 

In the following lemma we state some properties of the Bessel function (see [I] , [TU] for 

details). 

Lemma 2.1. The Bessel function J m (r) satisfies the following properties. 

J m (r) = 0(r m ), asr^O, (2.7) 
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N N 

J m (r) = e- ir ^a mj r- (j+ 5) + e ir ^a mij r- (i+ ^ +0(r- (iv+ f ) ), as r ->• +00, (2.8) 
3=0 3=0 

for any N G Z +; and 

n — 2 

r- — J^2(r) = c n Ji(e tr h(r)), (2.9) 
where h is a smooth function satisfying 

\d k r h(r)\<c k (l + r)-^- k , (2.10) 

for any integer k > 0. 

Finally, our last goal is to derive a Strichartz estimate. As usual, we will obtain 
it from a TT* argument, using an L 1 — L°° decay estimate for the linear propagator 
together with the conservation of its L 2 norm. 

We start by obtaining the decay estimate whose proof depends on the Van der Corput 
lemma. 

Lemma 2.2. Let f be a real valued C 2 function, and if) a C 1 function, defined in [a, b], 
such that |/" (01 > 1 f or an y £ £ [a,b]. Then 

\ txm m^\<\\\-H\m\ + u'\\L H1 a, b] )), (2.11) 

where the implicit constant is independent of a or b. 

The crucial decay estimate now follows. 
Lemma 2.3. Let < ft < 1, then for all t el and x G M 2 , one has 

\\DZu(t)w \\ LS > <r^|KIUi. (2.12) 

Proof. The idea of the proof of this lemma is very similar to the one used in [7j . Using 
the Fourier transform and Young's inequality we have 

mu{t)w \\ LT < \\(\ . iv^nMNiUi. (2.13) 
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The formula for the Fourier transform of a radial function (12.51) and the property (I2.9p 
yield 



-oo 

„3 



(|e|V^) v (x)= / s p e^ s+s \j (rs)sds 
Jo 

p+oo r+oo 

= / S P +1 e it{s3+s) e irs h(rs)ds + / s p+l e it{s " +s) e- irs Ws)ds 
Jo Jo 

= Ap(r,t) + Bp(r,t), 

(2.14) 

where r = \x\. Without loss of generality we consider t > 0. The estimates for both 
terms in (I2.14p follow the same techniques and the Ap(r,t) case is actually simpler, so 
we consider only the second term, Bp(r, t). Performing the change of variables u = t^s, 
we have 

/■-f-OO 2 1 

B p {r,t)=r*? / u p+1 e l ^ + ^ u - Tt ~* u) h(rt-\u)du 

Jo (2.15) 



where the function Mp{r) is defined as 

r+oo 

Mp(r) = / u p+1 e ifr Mh(ru)du, 
Jo 

2 

with the phase function, for 7 = £3, given by 

f r (u) = u 3 - (r - j)u. 
The proof of the lemma then reduces to proving that 



sup|M^(r)|<l, V/3e[0,l]. (2.16) 



r>0 

We now use standard oscillatory integral methods to obtain (I2.16p . By the stationary 
phase principle, the main contribution to the integral comes from the points where the 
derivative of the phase vanishes, i.e. 

3m 2 = r — 7. 

So we isolate the integral around these points, with a smooth cut-off function ip : M — > 
< if) < 1, rescaled around u = J ^-jp, such that 

2 /„ _m ^ r -7- 



supp^C{n:|3u - (r - 7)] < — — }, 
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and 

ip(u) = 1 for \3u — — 7)| < — - — . 

Notice, as this will be important below, that the length of the support is of the order 
(r — 7) 1 / 2 and that each derivative of ip produces a factor of (r — r ))~ 1 / 2 . Splitting the 
integral Mp(r) into two terms, with this cut-off function, we obtain 

\M p (r)\<\M}(r)\ + \M%(r)\, 

where 

r+00 

M 1 p{r)= U P +1 e iMu) h{ru)ifj{u)du, 
Jo 

and 

r+00 

Mj(r)= / u p+1 e ifr{u) h(ru)(l-ip(u))du. 
Jo 

When r — 7 < 0, no stationary phase points exist and no splitting is required, with 
the whole integral reducing to the M|(r) case. This is, in fact, what happens for the 
Ap(r,t) term above, whose phase is given by f r (u) = u 3 + (r + j)u and therefore has 
no stationary phase points in the domain of integration u G (0, 00). 

The bound for the non-stationary phase term Mj(r) is easily obtained integrating 
by parts in the usual way, using the non- vanishing derivative of the phase (see [7] for 
details). 

As for Mj (r), we now use Van der Corput's lemma to obtain an accurate bound 
for the integral around the stationary phase points. We start by writing f r (u) = 
(r— 7) 1 / 2 ^ ( r _" 3 )i/ 2 — (r— 7) 1//2 nj and note that on the support of ip we have u ~ (r— 7) 1 / 2 
so that 

v(r- 7 )V2 {r 1} a ) ( r _ 7 )i/2~ A » 

therefore 

H(r)\ < — ^(ll^rOCO^VCOlU- + \\^{Hr-)(-r +1 ^-))\\ L i). (2.17) 

Using the estimate for h and its derivatives (12.101) . the fact that on the support of ip 
we have u ~ (r — 7) 1 / 2 , and the observations made earlier about the derivatives of ij) 
and the length of its support, it is easy to check that all the terms on the RHS of the 
inequality (I2.17P are equally bounded by a constant multiple of 

( r _ 7 )(/3+l)/2 

(r- 7 ) 1 /4(i +r ( r _ 7 )i/2)i/2' 
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yielding 

H(r)\<l, 

with a constant depending on (3. □ 

Having obtained the decay estimate, we can now proceed to use it as one of the main 
ingredients in the proof of the Strichartz estimate for the unitary linear propagator. 

Theorem 2.3. Let T > and < 5 < §. Then, for all w G L 2 (M 2 ) 

\\DiUw \\ L « LS , < KHz*, (2.18) 

where q s = 

Proof. Proofs of these types of estimates, as consequence of the L 2 conservation and 
the L 1 — L°° decay estimates for the corresponding linear evolution groups, are now a 
standard procedure (see [2] or [3] for example). 
We start by interpolating 

\\U(t)w \\ L 2 = ||w ||l3, 

and 

\\D^U{t)w \\ Ls > <t-^\\w Q \\ L i, 

to obtain 

\\D^Mt>o\\ L i<t-^ 6 \\w \\ Lp/ , 

forO<e<l, I = I-f and 1 = J + ^. 

Then, using a TT* argument on the dual formulation and the Hardy-Littlewood- 
Sobolev inequality in the t variable, we get 

\\D^ 2 U(t)w \\ LtL , < \\w \\ L2 , 

for < 6 < 1, i = i - |, 2 - = ^6. 

— — ' p 2 2 ' q 3 

Finally, choosing the endpoint case 8 = 1 and writing 5 = (3/2 we obtain the desired 
inequality. Notice that, at 9 = 1, the case (3 = 1, i.e. 5 = \, becomes unattainable, due 
to the Hardy-Littlewood-Sobolev restriction ^-^-0 < 1. This corresponds to the known 
L 2 L^ Strichartz endpoint problem. □ 
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Before leaving this section we record two more results that will be required in the 
sequel. The first one is the Leibniz's rule for fractional derivatives whose proof can be 
found in [6]. 

Lemma 2.4. Let < m < 1 and 1 < p < oo. Then for any f, g : R" — > C, one has 

\\KU9)-fD™g-gD™f\\ Ll < \\g\\ L ~\\D™f\\ Ll . (2.19) 

The second result is a well known theorem in harmonic analysis, about the L p bound- 
edness of the Riesz transform, which states the following. 

Lemma 2.5. Let Ri, I = 1, 2 be the Riesz transform. Then, for 1 < p < oo 

WRifU < ll/IU- (2.20) 

3. Proof of the local well-posedness result 

In this section we provide a proof of the main result of this work. As explained 
in the introduction, we start by proving a local well-posedness result for the 4x4 
system that comes by differentiating each equation in (II. 3p with respect to X\ and re- 
introducing the new unknowns W\ := u Xl , W2 '■— u X2 , W3 := v Xl &nd := v Xa , we 
obtain the following IVP, 

'd tWl + zV=A(l - A) Wl = \(Rid x M + R 2 d Xl L 2 ) - \d Xl L 3 

d t w 2 + iV=E{\ - A)w 2 = \(Rxd x ,L x + R 2 d X2 L 2 ) - \d X2 L 3 

d t w 1 - ? V /Z A(1 - A)w 1 = ~l(R 1 d xl L l + R 2 d Xl L 2 ) - jd Xl L 3 

\ (3-1) 

d t w 2 - iyf^K(l - A)w 2 = ~l(R 1 d X2 L 1 + R 2 d X2 L 2 ) - \d X2 L 3 

w 1 (x 1 ,x 2 ,0) = wi(0), w 2 {x 1 ,x 2 ,0) = w 2 (0), 
w 3 (xi,x 2 , 0) = w 3 (0), w 4 (x 1 ,x 2 , 0) = w 4 (0) 

where 

Li := Lx{wi,w 2 ,w 3 ,w^) = (R1W1 + R 2 w 2 - Riw 3 - R 2 w i )(w 1 + w 3 ) 

L 2 := L 2 (wi 1 w 2l w 3l W4) = (R±wi + R 2 w 2 - Riw 3 - R 2 Wi)(w 2 + w A ) (3.2) 

L 3 := L 3 (w 1 ,w 2 , w 3 , Wi) = (w 1 + w 3 f + (w 2 + w 4 ) 2 . 

In what follows we give a local well-posedness result for the IVP (13. ip . Before 
presenting this result, we introduce some notations. Define J£ S (M 2 ) := (H S {R 2 )) with 
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norm given by ||fi'||^(R2) := E*=i lb.jlU*(R 2 ), for 9 = (9i, 92, 93, 9a) e JC S (M 2 ). We also 
denote by w = (wi,W2,Ws,w^) the space-time functions in (C([0,T] : fP(IR 2 ))) 4 = 
C([0,T] : IK S (M 2 )). 

Theorem 3.1. Let s > 3/2, then for any w(0) G IK S (IR 2 ), there exist a time T = 
T (\\w(0)\\3is(M. 2 )) > a space continuously embedded in C([0, T] : 3i" s (IR 2 )) and a unique 
solution w G to the Cauchy problem (13.1 p . 

Proof. We consider only the most difficult cases | < s < 2, as for s > 2 the proof is 
slightly simpler. Let the fractional part of s be m := s — 1 so that m G (1/2, 1), and 
let 5 = 1 - m = 2 - s, defining q s as in Theorem 1231 For / G C([0, T] : # S (M 2 )) define 



nUf) ■= sup H/IIh-, (3.3) 

0<t<T 

tf T (f) ■= E u#/ikzff + li^/iu^ + EE ll^/lk^ + E \\DlRifHL~, 

\P\<1 l=l,2\/3\<l 1=1,2 

(3.4) 



nUf) ■= E W&MWizlf + E E IPi^/ll^i*.. ( 3 - 5 ) 

|/3| = 1 i=l,2|/3| = l 



fi 4 (/) := E sup \\d%Dl + ™f\\ L% + E E SU P H W^lk x[0T1 > (3-6) 

n*(/) - ( E ll/llia.^)* + E ( E ( 3 - 7 ) 

«ez 2 J=l,2 aeZ 2 

and let ^lr(f) '■= m aXj=i,-,5 ^t(/)- 
Define the Banach space Xf. by 

X* := {/ G C([0,T] : H S (R 2 )) : fl T (/) < oo}, 

with norm ||/||xf, := ^t(/)j an d finally define := (Xf,) 4 with norm ||u>||xf, := 
E-=i IKIUj, identifying (C([0,T] : # S (R 2 ))) 4 with C([0,T] : H S (R 2 )). 
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For the given initial data u>(0) = (wi(0), w 2 (0), w 3 (0), w 4 (0)) G JC S (IR 2 ) and using 
Duhamel's formula, we define applications 



:= t/ + (t)(«;i(0)) + \ f* U + (t - ^{RAM + R 2 d Xl L 2 - \d Xl L 3 )(t')dt' 
q 2 ( W2 ) : = U+(t)(w 2 (p)) + | jl U + (t - t')(R 1 d X2 L 1 + R 2 d X2 L 2 - \d X2 L 3 ){t')dt> 
* 3 (^ 3 ) := U-(t)(w 3 (0)) - \ f*U~(t - ^{RxdvM + R 2 Xl L 2 + \d Xl L 3 )(t')dt 
My*) := U-(t)(w 4 (0)) - \ f*U~(t - ^{Rxd^Lx + R 2 d X2 L 2 + \d X2 L 3 )(t')dt> . 



Define a ball with center at the origin and radius a > in Xy by T> a := {/ G Xf* : 
ll/llxf, < a }- We will prove that there exist a > and T > such that the application 
\P := ^2, ^3, ^4) is a contraction map on the ball "B a . 

Using linear estimates (I2.2p . (I2.4j) . (j2.18j) . we have that 



*i( M i)IUf. < Co||tUj(0)|| fl «(R2) + C / HF^Wi,^, W3,W4)(t')||^(R2)^', (3.9) 



for all j — 1, • • • , 4 where Fj(w\, w 2 ,w 3 , W4) are the respective nonlinear parts in equa- 
tion f l3~Tj) . 

The estimates for each component ^(wj) follow similar techniques, therefore we 
present details only for the nonlinear part associated to the first component, ^i{w\). 
Note that in view of (13. ip and (13. 2p the terms in Fi(w) are of the form Rid Xl {(RiWj)wk} 
and d xl (wjWk), for j, k = 1, • • • , 4 and i, I = 1, 2. 

Again, we only provide details for the terms that are in the form Rid Xl {(RiWj)wk} , 
i.e., we find estimates for J Q T \\Rid xl {(RiWj)wk}(t')\\H=dt' . The estimates for the other 
terms follow a similar method and are even easier to handle. Using the fact that the 
Riesz transform is bounded in L p , 1 < p < 00, we have 



(3.8) 





(3.10) 



:/ + //. 



Now, 
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I< [ URid^w^^dt' + f \\{R lWj ){d Xl w k )\\ Ll dt' 
Jo Jo 

\\Rid xl Wj\\ L i \\w k \\ L ™dt' + / \\RiWj\\ L 2\\d Xl w k \\ L ™dt' 



< 



o JO 

3 

< \\w k \\L™L™T\\d xl Wj\\ L¥L 2 + \\wj\\L ¥ LiT^\\d xl w k \\ L 3 TLT 



(3.11) 



To estimate //, recall the notation a/— A = D\ = R\d Xl + R2d X2 so that we can 
write D s x = D\D™ = D™R x d Xl + D x a R 2 d X2 . Therefore we have 

With this decomposition, we can split II into four terms 

II< I \\D^(R l d Xl D 1 x w j w k )\\ L 2dt'+ [ ^(RA^DlwkjW^dt' 



JO 

T r T 



+ / \\D™(R l D 1 x w J d Xl w k )\\ Ll dt'+ [ WDZfawjDld^wJWvJlf (3 ' 12) 
Jo Jo 

=: a + b + c + d. 

Using Leibniz's rule for fractional derivatives ( I2.19p . L p -boundedness of the Riesz 
transform (12.201) and Holder's inequality, one gets 



< 



a ^ 



T W^Dlw^D^w^dt' + f UR^Dl^w^Wkhidt' 
o Jo 



< ||D> fc || L? , Li Tt||i? AlJ DX|l^^+^(E / / \w k Rid Xl Dl +m w 3 \ 2 dxdt) 



aez 2 



(3.13) 
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Likewise 

b< [ T \\d Xl R lWj \\ L ^\\Dl +m w k \\^dt' + f ' \\{D™d xl R lWj )Dlw k \\ Ll dt' 
Jo Jo 

^ \\D l x +m w k \\ L ^ L 2 \\d Xl Riw j\\ L oo dt' + WDlwkWL^WD^d^RiWjW^dt' ^.u) 

< Q^Wk^Wd^RiWjW^L^ + 1 1 D?d Xl | Ti 1 1 D x x w k \ \ L%L ~ 

< T^nKw^nKwj) + T^VL l T ( Wj )n 2 T (w k ). 

The estimates for the terms c and d are similar to those for the terms b and a, 
respectively, so we omit the details. Inserting the estimates fl3.10p - fl3.14p in f 13 . 9 1) . one 
obtains 

||*i WIU{. <co|K(0)||^ (R2) + C T £ (||H|xs) 2 , (3.15) 
for some e > 0. Gathering the similar estimates for the other components of we get 

II^HIIx^ <Co|k(0)|| MS(M2) +4cT £ (||«;|| x ^ 2 . (3.16) 



Given w(0) G !K S (IR 2 ), if we choose a = 2cQ\\w(0)\\wQjp) and T > in such a way 
that 4cT e ||w||xs i < ^, it is easy to show that $ maps 25 a onto itself. A totally analogous 
method is then used to prove that the application \1/ is a contraction mapping. The 
remaining part of the proof follows by a standard argument. □ 



Proof of Theorem li.il Having established Theorem 13.11 the proof of Theorem 11.11 fol- 
lows easily. For data i] G H S (M?) and <3> V S+1 (IR 2 ), we start by defining 

i 1 i 1 

U ° = ~27^E V ° + 2 $0 V ° = 2 v /z A r/ ° + 2 $ °' ^' 17 " > 

and then 

^i(O) = d Xl u , w 2 (0) = d X2 u , u;g(0) = ft^o, 104(0) = d^o, 

noting that, from our definition of V S+1 (R 2 ) we get w(0) G JC S (IR 2 ), with s > 3/2. 

We can now apply Theorem 13.11 to w(0), obtaining a unique solution w(t) G X^. 
Differentiating the first equation in (13 .ip in £2 and the second in x%, shows that 

d t (d X2 wi - d Xl w 2 ) = 0. 

Doing the same for the third and fourth equations, equally yields 

d t (d X2 w 3 - d xi w 4 ) = 0. 



BOUSINESQ TYPE SYSTEM 

From our construction of the initial data, we have 
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<9 X2 Wi(0) - d xi w 2 (0) = and d X2 w 3 (0) - d xl w A (0) = 0, 

which, from the two previous relations, proves that the vector fields (wi,W2) and 
(io 3 , W4} are closed in R 2 for all t G [0, T}. Therefore, there exist unique u, v G C([0, T] : 
V S+1 (IR 2 )) such that (d Xl u,d X2 u) = (wi(t),W2(t)) and (d Xl v, d X2 v ) = (w 3 (t), w 4 (t)), sat- 
isfying u(0) = Uq and v(0) = Vq according to (I3.17p . Finally, we invert back from (u, v) 
to (77, $) using formula ( II. 6p . Tracing the steps of this proof, we observe that this so- 
lution is obtained uniquely in the space ^ S T = {(77, $) : rj G Xf, \/ — A$ G Xf| which 
is continuously embedded in C([0,T] : ^(M 2 ) x V S+1 (M 2 )) = C([0,T] : ^(M 2 )) x 
C([0,T] : V S+1 (R 2 )). □ 



References 

[1] Z. Guo, L. Peng, B. Wang, Decay estimates for a class of wave equations, J. Funct. Anal., 254 
(2008) 1642-1660. 

[2] M. Keel, T. Tao, End-point Strichartz estimates, Amer. J. Math., 120 (1998) 955-980. 

[3] C. E. Kenig, G. Ponce and L. Vega, Oscillatory integrals and regularity of dispersive equations, 

Indiana Univ. Math. J., 40 (1991) 33-69. 
[4] C. E. Kenig, G. Ponce and L. Vega, Well-posedness of the initial-value problem for the Korteweg- 

de Vries equation, J. Amer. Math. Soc, 4 (1991) 323-347. 
[5] C. E. Kenig, G. Ponce and L. Vega, Small solutions to nonlinear Schrodinger equations, Ann. 

Inst. Henri Poincare, 10 (1993) 255-288. 
[6] C. E. Kenig, G. Ponce and L. Vega, Well-posedness and scattering results for the generalized 

Korteweg-de Vries equation via the contraction principle, Comm. Pure and Applied Math., 46 

(1993) 527-620. 

[7] F. Linares, D. Pilod, J-C Saut, Well-posedness of strongly dispersive two-dimensional surface 

waves Boussinesq systems, arXiv:1103.4159 ( 2011. 
[8] J. R. Quintero, Solitary water waves for a 2D Boussinesq type system, J. Partial Differ. Equ. 23 

(2010), 251-280. 

[9] J. R. Quintero, The Cauchy problem and stability of solitary waves for a 2D Boussinesq-KdV 
type system, Differential Integral Equations 24 (2011), 325-360. 
[10] E. M. Stein and G. Weiss, Introduction to Fourier Analysis in Euclidean Spaces, Princeton Uni- 
versity Press, Princeton, N.J., 1971. 



16 F. LINARES, M. PANTHEE, AND J. DRUMOND SILVA 

Felipe Linares: IMPA, Estrada Dona Castorina 110, 22460-320, Rio de Janeiro, RJ, 
Brazil. 

E-mail address: linares@impa.br 

Mahendra Panthee: CMAT, University of Minho, Campus de Gualtar, 4710-057, 
Braga Portugal. 

E-mail address: mpanthee@math.uminho.pt 

Jorge Drumond Silva: Center for Mathematical Analysis, Geometry and Dynamical 
Systems, Departamento de Matematica, Instituto Superior Tecnico, Av. Rovisco Pais, 
1049-001 Lisboa, Portugal. 

E-mail address: jsilva@math.ist.utl.pt 



